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ABSTRACT 

As a potential mechanism to build up supermassive black holes (BHs) in a spheroidal system, 
we consider the radiation drag effect by bulge stars, which extracts angular momentum from 
interstellar gas and thus allows the gas to accrete onto the galactic center. With incorporating 
radiation hydrodynamical equation with simple stellar evolution, it is shown that the BH-to- 
bulge mass ratio, /bhj is basically determined by a fundamental constant, that is, the energy 
conversion efficiency for nuclear fusion of hydrogen to helium, e = 0.007. More specifically, /bh 
is predicted to be 0.3e — 0.5e. Based on the present model for BH growth, a scenario for quasar 
formation is addressed in relation to ultraluminous infrared galaxies. 

Subject headings: black hole physics — accretion — galaxies: active — galaxies: formation — galaxies: 
evolution — galaxies: nuclei — galaxies: starburst 



1. Introduction 

The recent compilation of the kinematical data 
for galactic centers in both inactive and active 
galaxies has revealed that the estimated mass of 
a central 'massive dark object' (MDO), which is 
the nomenclature for a supermassive BH candi- 
date, does correlate with the properties of galactic 
bulges. The demography of MDOs has shown a 
number of intriguing relations: 1) The BH mass 
exhibits a linear relation to the bulge mass with 
the ratio of /bh = Af B H/Mbui ge = 0.001 - 0.006 
as a median value (Kormendy & Richstone 1995; 
Magorrian et al. 1998; Richstone et al. 1998; Geb- 
hardt et al. 2000; Ferrarese & Merritt 2000; Mer- 
ritt & Ferrarese 2001a) (It is noted that the bulge 
means a whole galaxy for an elliptical galaxy.) 2) 
The BH mass correlates with the velocity disper- 
sion of bulge stars with a power-law relation as 
M BH oc cr™, n = 3.75 (Gebhardt et al. 2000a) 
or 4.72 (Ferrarese & Merritt 2000; Merritt & Fer- 
rarese 2001a, 2001b). 3) The ratio / BH tends to 
grow with the age of youngest stars in a bulge until 
10 9 yr (Merrifield et al. 2000). 4) In disk galaxies, 



/bh is significantly smaller than 0.001 when the 
disk stars are included (Salucci et al. 2000; Sarzi 
et al. 2001). 5) For quasars, /bh is level with 
that for elliptical galaxies (Laor 1998). 6) As for 
Seyfert galaxies, the mass ratio is under contro- 
versy. The estimation by the reverberation map- 
ping suggests that /bh in Seyfert 1 galaxies is con- 
siderably smaller than 0.001 (Wandel 1999; Geb- 
hardt et al. 2000b), although the reverberation 
mapping method may cause a systematic under- 
estimate of BH mass (Krolik 2001). However, the 
measurements of Hp line widths give /bh — 0.0025 
for Seyfert l's as well as quasars (McLure & Dun- 
lop 2001). On the other hand, the BH mass-to- 
velocity dispersion relation in Seyfert 1 galaxies 
seems to hold good in a similar way to elliptical 
galaxies (Nelson 2000; Gebhardt et al. 2000b). All 
of these correlations imply that the formation of a 
supermassive BH is likely to be physically linked 
to the formation of a galactic bulge which harbors 
a supermassive BH. In addition, a great deal of re- 
cent efforts have revealed that quasar host galax- 
ies are mostly luminous and well evolved early- 
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type galaxies (McLeod & Rieke 1995; Bahcall et 
al. 1997; Hooper, Impey, & Foltz 1997; McLeod, 
Rieke, & Storric-Lombardi 1999; Brotherton et al. 
1999; Kirhakos et al. 1999; McLure et al. 1999; 
McLure, Dunlop, & Kukula 2000). These findings, 
combined with the BH-to-bulge relations, suggest 
that the formation of a supermassivc BH, a bulge, 
and a quasar is mutually related. 

Some theoretical models have been considered 
to explain the BH-to-bulge correlations, e.g., hy- 
drodynamical ones including a wind-regulation 
model (Silk & Rees 1998) or an inside-out ac- 
cretion model (Adams, Graff, & Richstone 2001), 
and self-interacting dark matter model (Ostriker 
2000). But, little has been elucidated regarding 
the physics on the angular momentum transfer 
which is inevitable for BH formation, since the 
rotation barrier by the tidal spin up in a growing 
density fluctuation is given by 

i^;~ 10 U^J (,005 ) (1+z) 

in units of the Schwarzshild radius i?s c h, where 
Mb is the baryonic mass, z is the cosmological 
redshift, and A is the spin parameter which pro- 
vides the ratio of circular velocity to velocity dis- 
persion of dark matter (Sasaki & Umemura 1996). 
Here, Rb arr is estimated by Rbarr = jb(z) 2 /GMf,, 
where jb is the specific angular momentum given 
by jb — -RmaxcA with i? max being the maxi- 
mum expansion radius of the fluctuation. Fur- 
thermore, required mechanisms for BH formation 
must work effectively in a spheroidal system like a 
bulge. The a- viscosity or non-axisymmetric gravi- 
tational instabilities would effectively transfer an- 
gular momentum once a disk-like system forms, 
but they are not likely to work in a spheroidal 
system. In this Letter, as a potential mechanism 
in a spheroidal system, the relativistic drag force 
by the radiation from bulge stars is considered, 
and the BH-to-bulge ratio is derived with incorpo- 
rating radiation hydrodynamics jointly with sim- 
ple stellar evolution in a bulge. As a result, the 
BH-to-bulge ratio is basically determined by the 
energy conversion efficiency for nuclear fusion of 
hydrogen to helium, e = 0.007. Also, in relation 
to BH growth, a scenario for quasar formation is 
addressed. 



2. Radiation Drag-Induced Mass Accre- 
tion 

We suppose a simple two-component system 
which consists of a spherical stellar bulge and 
dusty gas within it. The exact expressions for 
the radiation drag are found in some literatures 
(Umemura, Fukue, & Mineshige 1997; Fukuc, 
Umemura, & Mineshige 1997). For the total lumi- 
nosity i* of a uniform bulge, the radiation energy 
density is given by E ~ L*/cR 2 , where c is the 
light speed and R is the radius of the bulge. Then, 
the angular momentum loss rate by the radiation 
drag is given by dlnJ/dt ~ —\E/c, where J is the 
total angular momentum of gaseous component 
and x is the mass extinction coefficient which is 
given by \ = n/p with dust absorption coefficient 
k and gas density p. Therefore, in an optically- 
thin regime, dlnJ/dt ~ — Ti*/c 2 M s , where r is 
the total optical depth of the system and M g is 
the total mass of gas. In an optically-thick regime, 
the radiation drag efficiency is saturated due to 
the conservation of the photon number (Tsuribc 
& Umemura 1997). Thus, an expression of the 
angular momentum loss rate which includes both 
regimes is given by dlnJ/dt ~ — (L»/c 2 M g )(l — 
e~ T ). In practice, it is likely that optically-thin 
surface layers are stripped from optically-thick 
clumpy clouds by the radiation drag, and the 
stripped gas losing angular momentum accretes 
onto the center (Kawakatu & Umemura 2001). 
Since the radiative cooling is effective in the sur- 
face layers, the accretion is likely to proceed in an 
isothermal fashion until an optically-thick massive 
dark object forms. Then, the mass accretion rate 
is estimated to be 

M = -M s ^^ ( l-e-). (1) 

In an optically-thick regime, this gives sim- 
ply M = L*/c 2 , which is numerically M = 
0.1M Q yr _1 (L*/10 12 Lo). This rate is compara- 
ble to the Eddington mass accretion rate for 
a black hole with 10 8 M Q , that is, M^dd = 
O.2M yr^ 1 77" 1 (M B H/lO 8 M Q ), where ■q is the 
energy conversion efficiency. For the moment, 
j] = 0.42 for an extreme Kerr black hole is as- 
sumed. For the moment, an optically-thick stage 
is considered. Then, the mass of an MDO is esti- 
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mated by 

Mmdo = / Mdt ~ / L*/c 2 dt. (2) 
Jo Jo 

The next task is to construct a model for bulge 
evolution. Here, we employ a simplest analytic 
model. The gas fraction f g = M g /Mb of the sys- 
tem is regulated by 

%L = -S(t) + F{t), (3) 

where S(t) is the star formation rate which is as- 
sumed to be a Schmidt law, kf g , and F(t) is the 
mass loss rate from stars which includes super- 
nova explosions and stellar winds. If we invoke 
the instantaneous recycling approximation, then 
F(t) = S(t)(l — a), where a is the net efficiency 
of the conversion into stars after subtracting the 
mass loss. Then, we have an analytic solution 
as f g — e~ akt . Hence, the star formation rate 
is M*/Mj = S(t) = ke- akt . The radiation en- 
ergy emitted by a main sequence star over its life- 
time is 0.14e to the rest mass energy, where e is 
the energy conversion efficiency of nuclear fusion 
from hydrogen to helium, which is 0.007. Thus, 
the luminosity of the bulge is estimated to be 
L* = Q.lAeke~ akt M b c 2 . By substituting this in 
(2), we find 

Mmdo = O-Uea^M^l - G - akt ). (4) 

The term M(,(l — c~ akt ) represents just the stellar 
mass in the system, that is, Mb u i ge observationally. 
As a result, the MDO mass to bulge mass ratio is 
given by 

M M Do/M bu i ge = 0.3ea^5, (5) 

where ao.5 = a/0.5. It should be noted that the 
ratio is time- independent and also regardless of the 
star formation rate k. This result implies that the 
ratio is basically determined by the total number 
of emitted photons. 

3. BH Growth and Quasar Formation 

From an observational point of view, the mass 
of an MDO is often regarded as BH mass. But, 
in the present model they should be distinguished 
from each other if one strictly defines the BH 
mass by the mass inside the event horizon. In 



the present radiation drag model, it is not very 
likely that all the angular momentum of stripped 
gas is removed thoroughly. Hence, the MDO pre- 
sumably forms a massive and compact rotating 
disk due to a little residual angular momentum. 
Then, the BH is more likely to grow via accretion 
in the MDO, rather than the prompt collapse of 
the MDO. If the mass accretion within the MDO 
is driven by the viscosity, the AGN activity may 
be ignited. Supposing the mass accretion onto the 
BH horizon is limited by an order of the Eddington 
rate, the BH mass grows according to 

M BH = M e^/ tEdd , (6) 

where v is the ratio of the BH accretion rate to 
the Eddington rate, v — MBn/MEdd, and ^Edd 
is the Eddington time-scale, ^Edd = 1.9 x 10 8 yr. 
Here M is the mass of a seed BH, which could be 
a remnant of a massive population III star with 
1O 2 ~ 3 M (Carr, Bond, & Arnctt 1984; Nakamura 
& Umemura 2001), or an early formed massive 
BH with ~ 10 5 M Q (Umemura, Loeb, & Turner 
1993). When the BH growth is given by (6), it 
is delayed from the growth of the MDO given by 
(4). In Figure 1, this delay is schematically shown. 
The BH mass reaches Mmdo at a time t C ross when 
(6) equals (4). We can calculate t C ross to find 
Across ~ 10^ _1 £Edd ~ 10 9 ^ _1 yr, although it has 
weak dependence on k. As seen in Figure 1, dur- 
ing t < t CIOSS , the BH mass fraction /bh increases 
with time and therefore possibly increases with the 
metallicity of the system. At t > icross, almost all 
of the MDO matter has fallen onto the central 
BH, and therefore the BH mass fraction is settled 
to / BH = Mmdo /Mbuigo = 0.3ea ] = 0.002^. 

So far, the optically-thick regime has been as- 
sumed. However, before t cmss the interstellar gas 
may be blown out by a galactic wind and the sys- 
tem could be transparent. Recently, it is argued 
that the color-magnitude relation of bulges can 
be reproduced if a galactic wind sweeps away the 
gas at the epoch of a few 10 8 yr (Kodama & Ari- 
moto 1997). If the system becomes optically-thin 
at a wind epoch i w , the radiation drag-induced 
mass accretion practically stops owing to the re- 
duced efficiency of radiation drag due to the small 
optical depth. Then, the final MDO mass is 
settled to M MD o = 0.14£a- 1 M 6 (l - e - akt ™). If 
the timescale of star formation (ak)^ 1 is longer 
than t w , the final MDO mass can be reduced. 
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But, MMDo/-^buigo is time-independent as stated 
above. When all gas of the MDO eventually falls 
onto the BH, again /bh = 0.3eag.5 i s achieved at 

Across • 

Based on the present model for BH growth, a 
picture of quasar formation may be constructed 
as shown schematically in Figure 1. Recently, 
the Eddington ratio is estimated to be v « 0.1 
for quasars (McLeod, Ricke, & Storrie-Lombardi 
1999). Hence, t CIOSS « 10 10 yr and therefore 
Across S> i w . Since the bulge luminosity decreases 
as L* oc c~ akt , the stage at t < t w is a bright, 
optically-thick phase, which may correspond to 
a ultraluminous infrared galaxy (ULIRG). Then, 
a ULIRG (t < t w ) harbors a more or less active 
nucleus. This can be regarded as a model for a 
paradigm of the evolution of ULIRGs to quasars 
proposed by Sanders et al. (1998) and Norman 
& Scoville (1998). Even at t > t w , M BH still 
continues to grow until t cross , and therefore the 
AGN brightens with time if the Eddington ratio 
is constant. This optically-thin phase may cor- 
respond to quasar phenomena. Its early phase 
might be observed as a narrow line type 1 AGN 
(Mathur 2001). After the AGN luminosity (L AGN ) 
exhibits a peak at t CIOSS , it fades out abruptly, 
not only because the mass accretion rate drops 
promptly in the optically-thin phase, but also be- 
cause the energy conversion efficiency is in propor- 
tion to M/MEdd for an optically-thin advection- 
dominated accretion flow (e.g. Kato, Fukue, & Mi- 
neshige 1998). The later fading nucleus could be a 
low luminosity AGN (LLAGN) (e.g. Kawaguchi & 
Aoki 2001). As a consequence, /bh increases with 
increasing Lagn up to ,/bh = 0.3£C% 5. Then, / B h 
for quasars is predicted to be similar to that of 
ellipticals, although possibly a bit smaller. As for 
/bh for Seyfert 1 galaxies, naively it is expected to 
be level with that for quasars, but circumnuclear 
starbursts might predominantly regulate the mass 
accretion (Umemura, Fukue, & Mineshige 1998). 

Finally, we point out that further mass accre- 
tion could be induced also by the AGN luminos- 
ity. The mass accumulated by this self-induced 
accretion is maximally Mi„d = r/(l — T]) Mbh, 
if the successive induction is included. Hence, 
at later stages as t > t CIOSS , we have /bh = 
(Mmdo + Mi nd )/M bulgc . Thus, the BH-to-bulge 
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Fig. 1 . — A schematic sketch of black hole growth 
and a scenario for quasar formation. The abscissa 
is time and the ordinate is mass or luminosity in 
arbitrary units. Mb u i ge is the mass of stellar com- 
ponent in the bulge. Mmdo is the mass of the 
massive dark object (MDO). Mbh is the mass of 
the supermassive BH. and Lagn are the bulge 
luminosity and the black hole accretion luminos- 
ity, respectively. t w is the galactic wind timescale, 
and t cr oss is defined so that Mmdo = Mbh- 
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mass ratio could become 

/bh = O.l^a-^l - J?)" 1 = 0.5ea^l, (7) 

if an extreme Kerr BH [j] — 0.42) is assumed, 
whereas the self-induction factor (1 — f])^ 1 = 1.06 
for a Schwartzschild BH (ry = 0.057) brings no sub- 
stantial increase on /bh- In the long run, although 
leaving some uncertainties with evolution model 
including ao.5, /bh is predicted to be 0.3e — 0.5e. 

4. Discussion 

The radiation drag efficiency could be strongly 
subject to the effect of geometrical dilution 
(Umemura, Fukue, & Mineshige 1998; Ohsuga 
et al. 1999). If the system is spherical, the emit- 
ted photons are effectively consumed within the 
system, whereas a large fraction of photons can 
escape from a disk-like system and thus the drag 
efficiency tends to be considerably reduced. Al- 
though this is very qualitative, the geometrical 
dilution may be the reason why /bh is observed 
to be significantly smaller than 0.001 in disk galax- 
ies. For the same reason, /bh could be lower for 
a flattened bulge. 

Finally, the present mechanism may provide a 
physical basis for the state-of-the-art cosmologi- 
cal scenarios for the joint formation of galaxies 
and quasars (Hachnclt & Rccs 1993; Haiman & 
Loeb 1998; Kauffmann & Haehnelt 2000; Monaco, 
Salucci, & Danese 2000; Granato et al. 2001; 
Hosokawa et al. 2001). In cosmological studies, 
/bh has been hitherto parameterized so that it 
should match to the observed BH-to-bulge mass 
relation. Also, Monaco, Salucci, & Danese (2000) 
argue that, in order to account for the observed 
statistics of quasars and elliptical galaxies, there 
should be a time delay between the beginning of 
the star formation and the quasar bright phase. In 
Granato et al. (2001), the formation of early-type 
protogalaxies as well as quasars is successfully ac- 
counted for if it is assumed that the relationship 
between the BH mass and the host mass has been 
imprinted during the early phase of the quasar and 
host evolution. In the present paper, the BH-to- 
bulge mass relation is a function of time, but it has 
been found that in a quasar phase their assump- 
tion holds good. In the light of these points, the 
basic assumptions and results by previous cosmo- 
logical scenarios for quasar formation arc mostly 



consistent with the present model for quasar for- 
mation based on the BH growth. 

The author thanks A. Ferrara, N. Kawakatu, 
S. Mineshige, T. Nakamoto, K. Ohsuga, and H. 
Susa for fruitful discussions. Also, the author is 
indebted to the support by Center for Computa- 
tional Physics, University of Tsukuba. This work 
was supported in part by the Grant-in- Aid of the 
JSPS, 11640225. 



5 



REFERENCES 

Adams, F. C, Graff, D. S., & Richstone, D. 0. 

2001, ApJ, 551, L31 
Bahcall, J. N., ct al. 1997, ApJ, 479, 642 
Brothcrton, M. S., et al. 1999, ApJ, 520, L87 
Carr, B. J., Bond, J. R., & Arnett, W. D. 1984, 

ApJ, 277, 445 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 

Fukue, J., Umemura, M., & Mineshige, S. 1997, 

PASJ, 49, 673 
Gebhardt, K., et al. 2000a, ApJ, 539, L13 
— . 2000b, ApJ, 543, L5 

Granato, G. L., Silva, L., Monaco, P., Panuzzo, 
P., Salucci, P., Zotti, G. De, & Dansc, L. 2001, 
MNRAS, 324, 757 

Haehnelt, M. G., & Rees, M. J. 1993, MNRAS, 
263, 168 

Haiman, Z. & Locb, A. 1998, ApJ, 503, 505 
Hooper, E. J., Impcy, C. D., & Foltz, C. B. 1997, 

ApJ, 480, L95 
Hosokawa, T., Mineshige, S., Kawaguchi, T. 

Yoshikawa, K. & Umemura, M. 2001, PASJ, 

in press 

Kato, S., Fukue, J., & Mineshige, S. 1998, Black- 
Hole Accretion Disk (Kyoto: Kyoto Univ. 
Press) 

Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 
311, 576 

Kawaguchi, T., & Aoki, K. 2001, preprint 
Kawakatu, N., & Umemura, M. 2001, MNRAS, 
submitted. 

Kirhakos, S., Bahcall, J. N., Schneider, D. P., & 

Kristian, J. 1999, ApJ, 520, 67 
Kodama, T., & Arimoto, N. 1997, A&A, 320, 41 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 

581 

Krolik, J. H. 2001, ApJ, 551, 72 
Laor, A. 1998, ApJ, 505, L83 

— . 2001, ApJ, 553, 677 
Magorrian, J., et al. 1998, AJ, 115, 2285 
Mathur, S. 2001, AJ, in press (astro-ph/0107163) 
McLeod, K. K., & Rieke, G. H. 1995, ApJ, 454, 
L77 

McLeod, K. K., Rieke, G. H. & Storrie-Lombardi, 
L. J. 1999, ApJ, 511, L67 



McLure, R. J., Dunlop, J. S., & Kukula, M. J. 

2000, MNRAS, 318, 693 
McLure, R. J., Kukula, M. J., Dunlop, J. S., 

Baum, S. A., O'Dea, C. P., & Hughes, D. H. 

1999, MNRAS, 308, 377 
McLure, R. J., & Dunlop, J. S. 2001, MNRAS, in 

press (astro-ph/0009406) 
Merrifield, M. R., Forbes, Duncan A., & Terlevich, 

A. I. 2000, MNRAS, 313, L29 
Merritt, D., & Ferrarese, L. 2001a, MNRAS, 320, 

L30 

— . 2001b, ApJ, 547, 140 

Mineshige, S., Tsuribe, T., & Umemura, M. 1998, 
PASJ, 50, 233 

Monaco, P., Salucci, P., & Dancse, L. 2000, MN- 
RAS, 311, 279 

Nakamura, F., & Umemura, M. 2001, ApJ, 548, 
19 

Nelson, C. H. 2000, ApJ, 544, L91 
Norman, C. & Scoville, N. 1988, ApJ, 332, 124 
Ostriker, J. P. 2000, Phys. Rev. Lett., 84, 5258 
Ohsuga, K., & Umemura, M. 2001, ApJ, in press 

(astro-ph/0105474) 
Ohsuga, K., Umemura, M., Fukue, J., & Mi- 
neshige, S. 1999, PASJ, 51, 345 
Rees, M. 1984, ARA&A, 22, 471 
Richstone, D., et al. 1998, Nature, 395A, 14 
Salucci, P., et al. 2000, MNRAS, 317, 488 
Sarzi, M., et al. 2001, ApJ, 550, 65 
Sasaki, S. & Umemura, M. 1996, ApJ, 462, 104 
Silk, J., & Rees, M. 1998, A&A, 331, LI 
Tsuribe, T., & Umemura, M. 1997, ApJ, 486, 48 
Umemura, M., Fukue, J., & Mineshige, S. 1997, 
ApJ, 479, L97 

— . 1998, MNRAS, 299, 1123 
Umemura, M., Locb, A., & Turner, E. L. 1993, 

419, 459 
Wandel, A. 1999, ApJ, 519, L39 



6 



